home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SSIFA.z / SSIFA
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSSSSSIIIIFFFFAAAA((((3333FFFF))))                                                            SSSSSSSSIIIIFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SSIFA   - SSIFA factors a real symmetric matrix by elimination with
  10.      symmetric pivoting.
  11.  
  12.      To solve  A*X = B , follow SSIFA by SSISL.  To compute  INVERSE(A)*C ,
  13.      follow SSIFA by SSISL.  To compute  DETERMINANT(A) , follow SSIFA by
  14.      SSIDI.  To compute  INERTIA(A) , follow SSIFA by SSIDI.  To compute
  15.      INVERSE(A) , follow SSIFA by SSIDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE SSIFA(A,LDA,N,KPVT,INFO)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAA REAL(LDA,N)
  25.         the symmetric matrix to be factored.
  26.         Only the diagonal and upper triangle are used.
  27.  
  28.      LLLLDDDDAAAA INTEGER
  29.         the leading dimension of the array  A .
  30.  
  31.      NNNN INTEGER
  32.         the order of the matrix  A .  On Return
  33.  
  34.      AAAA a block diagonal matrix and the multipliers which
  35.         were used to obtain it.
  36.         The factorization can be written  A = U*D*TRANS(U)
  37.         where  U  is a product of permutation and unit
  38.         upper triangular matrices , TRANS(U) is the
  39.         transpose of  U , and  D  is block diagonal
  40.         with 1 by 1 and 2 by 2 blocks.
  41.  
  42.      KKKKPPPPVVVVTTTT INTEGER(N)
  43.         an integer vector of pivot indices.
  44.  
  45.      IIIINNNNFFFFOOOO INTEGER
  46.         = 0  normal value.
  47.         = K  if the K-th pivot block is singular.  This is
  48.         not an error condition for this subroutine,
  49.         but it does indicate that SSISL or SSIDI may
  50.         divide by zero if called.  LINPACK.  This version dated 08/14/78 .
  51.      James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab.  Subroutines and
  52.      Functions BLAS SAXPY,SSWAP,ISAMAX Fortran ABS,AMAX1,SQRT
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.